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TWO-STEP MECHANISM FOR 


ORDER-DISORDER KINETICS IN SILICATES 
ABSTRACT 

It is proposed that the apparent high activation energy barrier associated 
with the lower temperature limit T Xr of cation ordering in certain silicates is 
attributable to a mechanism involving two consecutive steps in the ordering 
process. This mechanism is regarded as responsible for certain ordering 
characteristics which distinguish metamorphic, igneous plutonic and volcanic 
pyroxenes. An important consequence is that order-disorder equilibrium can 
be attained at temperatures below T Tr only through recrystallization and not by 
thermal annealing which does not result in recrystallization. 
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TWO-STEP MECHANISM FOR 


ORDER-DISORDER KINETICS IN SILICATES 


INTRODUCTION 

The kinetics of ionic intracrystalline exchange (order-disorder phenomena) 
is an important aspect of the thermal behavior of minerals. The recent experi- 
ments of Virgo and Hafner (1969) have indicated that the high temperature 
kinetics of order-disorder in orthopyroxene are consistent with a simple exchange 
model which had earlier been found applicable to metallic alloys (Dienes, 1955) 
and which has been adapted to silicates (Mueller, 1967). However Virgo and 
Hafner (1969) also made the important observation that there is an apparent cut- 
off or transition region on the temperature scale below which no more ordering 
or disordering occurs. They attributed this behavior to an unmeasured energy 
barrier different from the one which is operative above the cut-off temperature 
of approximately 480 °C in orthopyroxene and which has a measured activation 
energy of about 20 kilocalories in the direction of disordering. 

It is apparent that the postulated energy barrier which is operative below 
the cut-off range is sufficiently high to prevent ordering even on a geological 
time scale since natural pyroxenes of slowly cooled metamorphic rocks still 
possess considerable disorder (Ghose, 1964). As was pointed out by Virgo and 
Hafner the 20 kilocalory barrier is not sufficiently high to prevent such ordering 
and cannot account for the relatively narrow cut-off temperature region. Further 
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since natural cummingtonites from metamorphic rocks also show considerable 
disorder (Ghose, 1961), it is likely that a similar cut-off or transition in 
mechanisms applies in this case and other cases also. It is obvious however 
that the transition temperature T Tr need not be the same in all silicates. 


MECHANISM OF CONSECUTIVE STEPS 

Two alternatives suggest themselves for the change of kinetic mechanisms 
which appears to be required to explain the ordering cut-off temperature T Tr . 
These are (1) a two consecutive step mechanism involving a low temperature 
process with a high activation energy and a high temperature process with a 
lower activation energy, and (2) two competing processes. We shall adopt the 
first alternative and show later why the second alternative does not apply. 

The exchange kinetics of species such as Mg 2+ and Fe 2+ in a quasi-binary 
crystal of the type of orthopyroxene may be represented by the following second 
order rate equation: 


dX t 

"dT 


|c 0 [k 12 0 12 x, (1 - x 2 ) - k 21 a 21 x 2 (1 - X,) 


(i) 


In this equation Xj and X 2 are the site occupancy fractions of either of the ionic 
species of a binary pair in the sites 1 and 2 respectively (represented in ortho- 
pyroxene by the symbols and M 2 ), K 12 and K 21 are the disordering and 
ordering rate constants, while 4> l2 and , are corresponding factors analogous 
to activity coefficient products in macroscopic systems, except that they also 
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contain the thermodynamic properties of the activated state. C 0 is the total 
concentration (in moles or atoms cm -3 ) of both sites 1 and 2. In the case of 
completely ideal solutions, including both the stable lattice sites and the activated 
state sites 4> 12 ~ $ 2 i = ■*- Eq. (1) reduces to its simplest form. 

The potential energy relations which are reflected in the experimental 
activation energies and which are the chief contributors to the temperature 
sensitivity of the K 4 } factors may be represented as in Fig. 1. This figure in 
particular is intended to refer to the exchange process above the temperature 
T Xr as manifested in such experiments as those of Virgo and Hafner (1969). hi 
the latter case the experimental activation energies E * x 1 2 and E e * x21 which 
correspond to AU * 2 and AU* x for disordering and ordering were found to be 

E e * xl2 = 20 kilocalories 

E ex2 i - 16 kilocalories 

The cut-off temperature T Tr can be best explained by resorting to two 
consecutive processes with potential energy configurations as depicted in Fig. 2. 
In this figure the high temperature (H) activated state corresponds to the 
20 kilocalory barrier discovered in orthopyroxene by Virgo and Hafner while 
the low temperature (L) activated state corresponds to the barrier assumed to 
be operative below the cut-off temperature T Xr . According to kinetic theory 
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(Glasstone, Laidler and Eyring, 1941} the system must pass over the highest 


barrier but need not pass through the intermediate state. 

The mechanism of two consecutive reactions cannot however be explained 
solely on the basis of potential energy diagrams such as Fig. 2 or the corre- 
sponding experimental activation energies. Kinetic theory tells us that the 
absolute reaction rate constants K 12 and K 21 are determined by the Gibbs free 
energy of activation AG* (Glasstone, Laidler and Eyring, 1941) since 

k T 

K = v — r~ exp (- AG*/RT) . (2) 

In this expression v is the transmission coefficient, k is Boltzmann's constant, 
h is Planck's constant, and R is the gas constant. An additional relation is 

AG* = ah* - TAS* , (3) 

where AH* and As* are the enthalpy and entropy of activation respectively. In 
condensed systems such as this AG* will not differ greatly from AF* the Helm- 
holtz free energy of activation and AH* is consequently very similar to AE* the 
internal energy of activation. Also these latter values should be similar in 
magnitude to the experimental activation energy E* x . 

The transmission coefficient in Eq. (2) is difficult to evaluate and in the 
absence of any evidence to the contrary is usually set equal to unity. However 
it is clear from Eqs. (2) and (3) that even if AH*, AE* and E e * x are small AG* 
may be indeterminately large depending on T and AS*, assuming that AS* < 0. 
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The possible effect of different values of AG* and E g * x on the rate constants 
and their temperature sensitivities is depicted in Fig. 3. Since the initial state 
is the same for both high and low temperature processes and since the high 
temperature process appears to be governed by the simple second order rate 
Eq. (1), we may assume that the equation for the low temperature process has 
the same form and that only the values of the K$ factors of the two processes 
will differ throughout most of the temperature range. However it is postulated 
here that the K<fi factors for the two processes have the same value at the 
cut-off temperature T t r . 

hi consecutive processes the over-all rate is governed by the step for which 
K fp has the lowest value. Consequently below T Tr , K L < K H c/- H while above 
this temperature K n < K L <£ L . We may also imagine another diagram analogous 
to Fig. 2 but with AG* plotted instead of the potential energy. In such a diagram 
the height of the barrier would change drastically with the temperature so that 
in the high temperature range the H barrier would be higher than the L barrier. 
In this way it is possible to explain why at these temperatures only parameters 
such as E e *” which refer to the high temperature process are registered. It is 
interesting that the negative value of the entropy of activation which is required 
for the high temperature process has been postulated on other grounds (Mueller, 
1969). 
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Below the cut-off temperature where K L 0 L is in control the temperature 
sensitivity of this rate constant is so great that it decreases rapidly to such 
low values that the exchange process is frozen even on the cooling time scale 
of metamorphic rocks. Thus it is possible that in a temperature interval of less 
than 10G°K the value of K, decreases some orders of magnitude. Although no 
data are as yet available for a quantitative calculation, it is possible to get a 
qualitative impression of the effect by substituting certain values of K L and 
the temperature interval. Thus from the general relation 


BhK L 

hT 


E 


*L 
e x 


KT 2 


(4) 


where E is regarded as constant we obtain the integrated expression 



(5) 


where = K R at T" - T Tr . 

If for any temperature interval AT = T" - T' = 100 °K and T" = T Tr = 

753 °K (480°C), and if we take (K'^/K^) = 10 6 we find that E e * x L % 140,000 calories. 
It is obvious that in a smaller temperature interval or for a larger change in K L 
E g * x j would be even larger. 

EXCLUSION OF THE COMPETING MECHANISM 

In a mechanism involving competing processes the governing process of 

the over-all rate is that which occurs most rapidly. This situation is depicted 
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in Fig. 4, vhich is a schematic drawing analogous to Fig. 3. It is clear from 
this figure however that a mechanism involving competing processes is incapable 
of explaining the observations in any simple way since the low temperature 
process involves a less steep slope below T Tr than that of the high temperature 
process. Consequently it is impossible for a great decrease in ordering (and 
disordering) rates to occur below T Xr . 

NATURE OF THE LOW TEMPERATURE PROCESS 

So far we have merely interpreted the observed ordering cut-off temperature 
as resulting from a two step mechanism without identifying the atomistic be- 
havior responsible. Because of the inferred high activation energy for the low 
temperature process it seems likely that bond breaking is involved. Perhaps 
also this process involves the activation of vacancies which are inferred to play 
aii important role in particle migration in crystal lattices (Dienes and Damask, 
1965). Unfortunately it is impossible to characterize the kinetics in any greater 
detail than this now and it is likely that this will be possible only after extensive 
experimental work hi the temperature range of T Tr . 

GEOLOGICAL CONSEQUENCES 

Some consequences of the experimentally determined kinetic behavior of 
orthopyroxene have already been pointed out by Virgo and Hafner (1969), but 
others come to our attention by considering the theoretical model. The general 
behavior of such silicates may best be summarized by a schematic diagram 
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such as that of Fig. 5 . This figure illustrates the extent of attainment of hetero- 
geneous and homogeneous equilibrium in the variety of environments encountered 
in natural samples. 

It is clear first of all from Fig. 5 that the macroscopic intercrystalline 
exchange equilibrium (line a) may in principle be attained and quenched in at any 
temperature at which crystallization or recrystallization can occur and merely 
reflects the latest of such events. However the establishment of this equilibrium 
requires the extensive reconstruction of the crystals involved except when un- 
mixing occurs on dimension scales of some tens of microns as a maximum 
(Mueller, 1969). For example an intercrystalline distribution established at 
solidus temperatures in a plutonic igneous rock can only be destroyed by re- 
crystallization at a lower temperature such as might be induced by shearing 
stresses. On the other hand the corresponding- intracrystalline exchange 
equilibria of such rocks (line b) will not ordinarily be quenched in at any 
temperature because of the slow cooling rate. The same conclusion also applies 
to metamorphie rocks (line c) above the temperature T Xr . However in these 
rocks the attainment of the intracrystalline equilibria is possible below T Tr if 
crystallization or recrystallization occurs at these temperatures. The basis for 
this conclusion is that the establishment of the intracrystalline distributions 
should generally occur more readily than those between crystals (Mueller, 1969). 

There is of course a slight probability that some igneous plutonic and 
metamorphie silicates could have their intracrystalline distributions quenched in 
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above the temperature T Tr , as is indicated by the dashed lines. This could 
occur for example if a recrystallizing metamorphic rock were cooled rapidly 
by introduced water along a fracture. 

Only in volcanic rocks however is there a high probability for quenching of 
the intracrystalline distributions and such distributions were indeed measured 
by Ghose and Hafner (1967). This characteristic is indicated by line d in 
Fig. 5. 

The lower segment of line c is particularly interesting because it indicates 
that if a slowly cooled crystal such as a metamorphic orthopyroxene crystallized 
or recrystallized above T Xr the observed distribution of Mg 2+ and Fe 2+ between 
the M sites all approximate the equilibrium distribution at T Xr except that T Xr 
itself is likely to vary somewhat as a function of compositional details. Thus 
it is unlikely that a suite of metamorphic orthopyroxenes define an equilibrium 
isotherm in detail. 
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Figure 1. Potential energy diagrams for an intracrystal I ine exchange process involving species 
A and B between sites 1 and 2. The shapes of the wells and barrier have no significance. 
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Figure 2. Potential energy diagram for the inferred consecutive two-step exchange process. The L-energy barrier is assumed to be 
in control at low temperatures while the H-barrier is assumed to be in control at high temperatures. The latter barrier also corresponds 
to that of Fig. 1. 
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Figure 3. Suggested relation between the rate factors K <fi of the high and low temperature steps 
of the two consecutive step mechanism and the temperature. The full lines indicate the tempera- 
ture ranges in which the rate factors govern the over-all rate. 
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Figure 5. Scheme illustrating the possibility of crystals attaining and retaining equilibrium dis- 
tributions. The full lines indicate temperature ranges for different materials in which such dis- 
tributions are readily attained and retained. Dashed lines indicate temperature ranges in which 
these possibilities are greatly diminished. 
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